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Abstract 

Local boundary conditions for spinor fields are expressed in terms of a 1- 
parameter family of boundary operators, and find applications ranging from 
(supersymmetric) quantum cosmology to the bag model in quantum chro- 
modynamics. The present paper proves that, for massless spinor fields on 
the Euclidean ball in dimensions d = 2,4,6, the resulting (^(0) value is in- 
dependent of such a 6 parameter, while the various heat-kernel coefficients 
exhibit a ^-dependence which is eventually expressed in a simple way through 
hyperbolic functions and their integer powers. 
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I. INTRODUCTION 



The choice of boundary conditions in the theories of fundamental interactions has always 
attracted the interest of theoretical physicists, not only as a part of the general programme 
aimed at deriving the basic equations of physics from a few guiding principles [1-8] , but also 
as a tool for studying concrete problems in quantum field theory and global analysis [9-11]. 

In particular, we are here interested in studying local boundary conditions for massless 
spin-| fields, whose main motivations may be summarized as follows [12,13]. 

(i) The Breitenlohner-Freedman-Hawking [14,15] boundary conditions for gauged super- 
gravity theories in anti-de Sitter space are local and are expressed, for spin-| fields, in 
terms of a projection operator. The rigid supersymmetry transformations between massless 
linearized fields of different spins map classical solutions of the linearized field equations, 
subject to such boundary conditions at infinity, to classical solutions for an adjacent spin, 
subject to the same family of boundary conditions at infinity [12]. 

(ii) In simple supergravity the spatial tetrad and a projection formed from the spatial com- 
ponents of the spin- 1 potential transform into each other under half of the local supersym- 
metry transformations at the boundary [16]. The supergravity action functional can also 
be made invariant under this class of local supersymmetry transformations. On considering 
the extension to supergravity models based on the group 0{N), the supersymmetry trans- 
formation laws show that, for spin-^ fields only, the same projector should be specified on 
the boundary as in the Breitenlohner-Preedman-Hawking case. 

(iii) The work in Ref. [13] has shown that, instead of quantizing gauge theories on a sphere 
or on a torus one can quantize them in an even-dimensional Euclidean bag and impose 
S'[/^(iVj)-breaking boundary conditions to trigger a chiral symmetry breaking. On investi- 
gating how the various correlators depend on the parameter 9 characterizing the boundary 
conditions one then finds that bag boundary conditions are a substitute for small quark 
masses [13]. 
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More precisely, in theories of Euclidean bags, chiral symmetry breaking is triggered by 



imposing the boundary conditions ||T3|,|T7| 

= vr.^loA/ = ^ (l + 



dM 



on the spinor field ip. Here we focus on the d- dimensional Euclidean ball, which is the portion 
of flat dimensional Euclidean space bounded by the S'^~^ sphere. The eigenspinors of the 



Dirac operator on the ball have the form ||18 

tJn+d/2{kr)Z^^\n) 



c 



^{d-2)/2 



.eJnHd-2)/2{kr)Zi\^) 
^Jn4-d/2(A;r)Zf)(^]) 



(J / eJn+(d-2)/2{kr) Zi"^ (fi) 



;i.2) 



:i.3) 



where C is a normalization constant, e = ±1, n = 0, 1, 2, oo, and Z^^\^l) are the spinor 
modes on the sphere |T9|. In Eq. (1.1), the boundary operator reduces to the matrix 

1 / 1 -te' 

2 i te-' 1 



and its application to (1.2) and (1.3) yields the eigenvalue condition ||20 



Jn+d/2{k) - ee^ Jn+d/2~l{k) = 



:i.4) 



for V'±'\ and 



Jn+d/2{k) +ee Jn+d/2-l{k) = 



:i.5) 



for ip^. ■* , where r has been set to 1 for convenience [3] . By eigenvalue condition we mean the 
equation obeyed by the eigenvalues by virtue of the boundary conditions, which yields them 
only implicitly [3]. Equations (1.4) and (1.5) lead eventually to the eigenvalue condition in 
non-linear form, i.e. 



Jn+d/2-lik) e '^^ Jn+d/2{k) — 0, 



:i.6) 



J:. 



n+d/2 



_i(/.)-e^Vt,/2(A:) = 0. 



;i-7) 
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Of course, it is enough to deal with one of these equations, while the contributions from the 
other follow by replacing 9 with —9. 

Recently, in ci = 2 dimensions, the spectral asymmetry following from the boundary con- 
ditions (1.1) was considered in [^. Asymmetry properties are encoded in the eta function 
which was analyzed using contour integral methods, see e.g. [TI[ . 



Instead, we study heat-kernel asymptotics for the squared Dirac operator on the d- 
ball with eigenvalue conditions (1.6) and (1.7) which is related to an analysis of the zeta 
function. Strictly, one can actually obtain two second-order operators of Laplace type out 
of the Euclidean Dirac operator i.e. 

Pi = DD^ and Pa = D^D, 

where P*^ denotes the (formal) adjoint of D. The existence of both Pi and P2 is crucial 
for index theory [5], and by taking into account both (1.6) and (1.7) we correctly take care 
of this (see Ref. [12] for the mode-by-mode version of Pi and P2 on the 4-ball). To be 
self-contained, recall that, given the second-order elliptic operator P, the heat kernel can be 
defined as the solution, for r > 0, of the associated heat equation 

(J:^ + P^Uix,y;T) = 0, (1.8) 

subject to the initial condition {{M,g) being the background geometry) 

lim / U{x,y;T)ip{y)Jdetg dy = (fix), (1.9) 

and to suitable boundary conditions 

[BUix,y;T)]9M = 0, (1.10) 

which preserve ellipticity and lead to self-adjointness of the boundary- value problem [9-11]. 
The functional (or L^) trace of the heat kernel is obtained by considering the heat-kernel 
diagonal U{x,x;t), taking its fibre trace TtvU{x,x;t) (since U{x,y;T) carries (implicit) 
group indices in the case of gauge theories), and integrating such a fibre trace over M, i.e. 



TTL2e-^^ = f TYvU{x,x]T)Jdet g dx. (1.11) 
Jm 

The asymptotic expansion we are interested in holds for t —>■ 0^ and has the form [9] 

oo 

Tr^.(e--^) ~ T-'/' ^"/'«n/2(P,B), (1.12) 

n=0 

where the heat-kernel coefficients a„/2(-P, -B) are said to describe the global (integrated) 
asymptotics and consist of an interior part c„/2(-P) and a boundary part &„/2(-P, B), i.e. 

a,/2(P, B) = Cn/2{P) + &n/2(P, B) . (1.13) 

At a deeper level, we might introduce a smearing function and consider instead the 
trace of fe~'^^, with / a smooth function on M. This takes into account the distributional 
behaviour of the heat kernel from the point of view of invariance theory (here "invariance" 
refers to the invariants of the orthogonal group, which determine completely the functional 
form of an/2{P,B) [9]). However, mode-by-mode calculations like the ones we are going to 
consider can be performed without exploiting the introduction of /, and hence we limit 
ourselves to using Eqs. (1.11)-(1.13). Section II describes the (^-function algorithm in Ref. 
2^ on the Euclidean (i-ball [Q, and Sec. Ill generalizes the work in Ref. [12] by showing 



that, on the 4-ball, non-vanishing values of 6 in Eqs. (1-6) and (1.7) do not affect the 
conformal anomaly. The hardest part of our analysis is then presented in Sees. IV and 
V, where heat-kernel coefficients are studied for arbitrary dimension d, with several explicit 
formulae in d = 2,4,6. Concluding remarks and open problems are described in Sec. VI, 
while relevant details can be found in the Appendix. 

II. THE MOSS ALGORITHM FOR THE L>-BALL 

The starting point in our investigation of the eigenvalue condition (1.7) for the purpose 
of heat-kernel asymptotics is the use of the ^-function at large x, which was first described 
in Ref. with application to 4-dimensional background geometries. However, since we 



are interested in the Euclidean rf-ball, we put no restriction on the dimension of M, denoted 



by d as in Sec. I, and we follow the general procedure as outlined by Dowker p3|. First 
we point out that on replacing the eigenvalues A„ of P by A„ + {x being a large real 
parameter), one has the ^-function at large x in the form 

({s, x') = 5:(A„ + x')-^ = —- / T^-'U^{r)dT, (2.1) 

n r(s) Jo 

having defined the integrated heat kernel (or functional trace of the heat kernel at large x) 



as 



f/,(r) = J2 e-(^"+"')" = e-"'"?7(r). (2.2) 

n 

By virtue of the asymptotic expansion already encountered in the Introduction, i.e. 



U{t) = J2 e-'"^ ~ E «n/2r^ as r ^ 0+, (2.3) 

n n=0 



we therefore find 



CO 

C(s,a;^) ~ — - ^ a„/2/(x; s,n,ci), (2.4) 

[^) n=0 



having defined 



I(x;s,n,d)= / r^-^+^e"" " dr. (2.5) 
Jo 

Now we distinguish two cases, depending on whether d is even or odd. In the former, we 
consider s = s such that s — 1 — | = 0, i.e. s = 1 + | which implies (on defining rx^ = z) 

I{x- s, n, d) = rte""'" dr = x'^-^r (l + ^) , (2.6) 

and hence yields, for d = 2k, k = 0,1, 2, ... 

In the latter, we consider s = s such that 

s = l+ ^ . ^ = s, 



which imphes 

I{x; s, n, d) = r'^e-^'^ dr = x'^-^r (^^) ■ (2.8) 

On the other hand, since the function expressing the eigenvalue condition (1.7) admits 
a canonical product representation (see Appendix), one can prove, on setting i/ = n + | for 
d even, the identity 

r(..C<(..J..) = ,-.,.|-..-f"-)(i£)"' 

X \og[{ix)-'^^-'\jl_,{ix) - e^'Jliix))], (2.9) 

where 22 is the dimension dg of spinor space, and deg(n) = \ds I is the degener- 

\ n I 

acy associated with the imphcit eigenvalue condition (1.7). Thus, the heat-kernel coefficient 
ai is equal to ^^^^ , ^ (respectively ^p+ry) times the coefficient of x"'"^ (respectively a;~^~^) 
in the asymptotic expansion of the right-hand side of (2.9) in even (respectively odd) di- 
mension. On focusing for definiteness on the even d case, we now exploit the identity 

J,^i(k)^Jl(k) + '^Mk), (2.10) 

and obtain 



2u 

,x^ I ' ix 



Jl_Sx) - e^'^Jl{ix) = Jl^{ix) - I ^ + e'^ ) Jl{ix) + ^J,(ix)Jl(ix). (2.11) 



Thus, on defining ai, = y/u^ + x^ and using the uniform asymptotic expansions of Ji,{ix) 
and Jl{ix) summarized in the Appendix we find 

Jl_S^) - e^'^Jliix) ~ M!il^a,e2'^-e-2'^'°s{-+"'') \^lAe{t) + + 2iEiE2l , (2.12) 

27r L J 

where we have defined 

t=—, (2.13) 



Ae{t) = l + {e-l){l-e^'). 
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(2.14) 



As expected, our formulae reduce, at ^ = 0, to the asymptotic expansions used in Ref. [12]. 
From now on we need to recall that the functions Si and S2 have asymptotic series in the 
form 

Si-E^, (2.15) 

k=o 



S.~E#, (2.16) 



where Uk and are the Debye polynomials given in Ref. The asymptotic expansions 



on the right-hand sides of (2.15) and (2.16) can be re-expressed as 

Ukit) ~ a jit) 



E^-Er^, (2.17) 



k=0 



=0 M 



f^l.fMI. (2.18) 



where 



a.(t) = 6,(t) = > 0. (2.19) 

Now the asymptotic expansion (2.12) suggests defining 

S = Y.\Ae{t) + E2 + 2tSiE2, (2.20) 

and hence studying the asymptotic expansion of log(S) in the relation to be used in (2.9), 
i.e. 

~ - log(27r) + log + 2a^ - 1v log(z/ + a^) + log S. (2.21) 
^From the relations (2.13)-(2.20) S has the asymptotic expansion 



00 



S ~ E (2.22) 



where the first few Cp coefficients read 



Co^Ag + l + 2t, 

ci = 2aiAo + 26i + 2t(ai + 61), 
C2 = (2a2 + al)Ae + (262 + bj) + 2t{a2 + 62 + Qi&i), 

C3 = 2(a3 + 0102)^0 + 2(63 + 6162) + 2t{as + h + 0162 + a2&i). 
Now defining 

Co 



and making the usual expansion 



00 rj 

iog(i+/)=E(-iy^ T' 

j=o ■' 



vahd as / — > 0, we find 



logE = log Co + logE ~ log Co + E — 

p=i I'^j'J 

where explicit formulae can be given for all Ap. In particular 

Co 



Co 2 



A3 = ^-AiA2-i(Ai)3. 

Co 6 



Using the definition 



jointly with our previous formulae and the explicit form of Debye polynomials pj] , a lengthy 
calculation yields as many Ap terms as are needed. For example (cf. Ref. [12]) 



where 



fio(t) = l-lt' + lt\ (2.37) 



^iW = -^ + |^' + ?-fi', (2.38) 



1^2^ = -- + -t' - -t\ (2.39) 

25 531 2 221 4 1105 ^ 

cuo t = + 2.40 

^ ^ 512 320 64 576 ^ ^ 

^ 16 16 4 8 16 16 8 ' ^ ' 

uj2{t) = -L + Le + ---t^-- + -f, (2.42) 

^ ^ 16 16 8 8 16 16 ^ ^ 

1 ^4 .6 

^^^t) = -- + --- + -. (2.43) 

^ 24 8 8 24 ^ ^ 

In the calculation, all factors 1 + t in the denominators of Ai,A2 and ^3 have cancelled 
against factors in the numerators, as in the 6 = case [12]. Moreover, a simple but non- 
trivial consistency check shows that, at ^ = 0, equations (2.34)-(2.43) yield Ax^A^ and A3 
in agreement with Ref. [12]. 
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As a result of all these formulae we find 

oo 

\og[{ix)-^^''-^\jl_, - e^'Jl){ix)\ ~ 5:^.(i.,«.(a:)), (2.44) 

i=l 

where the first few functions Si read (cf. Ref. [12]) 

^1 = -log7r + 2a„ (2.45) 

^2 = -{2v - 1) log(z/ + a,) + log/e(t), (2.46) 
^3 = — , (2.47) 
^4 = ^, (2.48) 

S, ^ ^. (2.49) 

Expanding deg(?T,) in powers of ra, the infinite sum over n in the expression (2.9) can be 
evaluated with the help of formulae derived using contour integration, i.e. p2| : 



p=0 



2r (A: + I) 



(2.50) 



2r(i + i) ~'2- 

Here / is a real number larger than 1 and Bq = 1,-82 = |j-B4 = — etc. are Bernoulli 
numbers. In arbitrary dimension d, the expansion of deg(r;,) in powers of n is cumbersome 
and a systematic formula suitable for all d is given in Eq. 



III. CONFORMAL ANOMALY ON THE 4-BALL 



As a first application, we show how the ({0) calculation of Ref. [12] is extended to 
our boundary conditions involving 6, leading in turn to the eigenvalue conditions (1.6) and 
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(1.7). Only calculations for (1.7) are presented, but the full C(0) value, which expresses 
the conformal anomaly for a massless Dirac spin-^ field (we do not study the coupling of 
spinor fields to gauge fields, which would lead us instead to the subject of chiral anomalies), 
receives a contribution from (1.6) obtained by replacing 6 with —6 in the result from (1.7). 

In 4 dimensions. Sec. II shows that C(0) = 02 is equal to \ times the coefficient of x"^ in 
the asymptotic expansion of the right-hand side of Eq. (2.9) at d = 4. On setting m = n + 2 
the latter reads (here am{x) = y/m^ + x^) 

00 / ^ rl \ ^ °° 

^(m^-m) -- Y.Sr{m,<^m{x))r.Y.Wl, (3.1) 
m=o \zxaxj .^^ j^-^ 

with corresponding to the third derivative of S^, for all i. The terms contribute to 
02 in 4 dimensions only up to i = 5, and hence only their analysis is presented hereafter. 



2 

00 



Contribution of and 

The term is given by 

00 / 1 r/ \ ^ 

= E - ^) (^^j (-log(^) + 2«^)' (3-2) 

which is unaffected by ^-dependent boundary conditions. Thus, we know from Ref. [12] 
that Eqs. (2.50) and (2.51) imply vanishing contribution to C(0)- 
The term reads (here t = —) 

\' 

W^^ + W'^^, (3.3) 



00 (id 

wl=Y. - ^) Tr:r 1 - (2^ - 1) iog(m + «„) + log fe{t) 



where W"^^ is the 6'-independent part while denotes the part involving log fe{t). ^From 
Ref. [12] we know that contributes 

(M(o) = L + J_ ^ J_. (3.4) 

^ ^ ^ 120 24 30 ^ ' 

The log fe is dealt with by defining 

7 = ^(l-e''^) = -e'^sinh^, (3.5) 
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and hence writing (see (2.33)) 

log/, = log(l - 7) + log(l + /^ma-i). (3.7) 

At this stage we can exploit Eq. (2.28), with / = (3ma^ < 1 since is always evaluated 

at large x, and hence we find 

1 00 00 
W'J = ^ E(-l)'(^ + 2)(^ + 4)/3'= E (^'^' - m'^'h-'-', (3.8) 

° k=l m=0 

where the interchange of the orders of summation is made possible by uniform convergence. 
Interestingly, this sum contributes infinitely many terms with equal magnitude and 
opposite sign, so that ^(0) is unaffected. More precisely, we consider odd and even values of 
k and hence define 

00 00 
^ E ^2 = E m''^'o.-''-\ k^O,l, 2, (3.9) 

m=0 m=0 

Fs^± m^'^'a-^''-\ F4 ^ E m''^'o.-n^'-\ k ^ 1,2, ... (3.10) 

m=0 m=0 

By virtue of (2.51), Fi contributes to with zero weight for all k because of the cos ^ 
coefficient. Moreover, F2 is proportional to x~^ by virtue of (2.50). The sum F3 is instead 
proportional to x"^ (again by (2.50)), while F4 is such that its contribution Sp^ix] k) to C(0) 
reads 

2) - h.ix; 1) = (-IB + ^^H) = 0, (3.12) 

and infinitely many other relations along the same lines. 

Effect of Wl,Wt^ andWl 

The term is equal to 
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\ m "1 / 

= ly^^ + ly^^ + w^^^ + w^^^, (3.13) 

where VF^"^ and are the first two, ^-independent sums, while W^f^ and are the 

sums depending on Q through 7 and /3. By virtue of the large-x nature of the whole analysis, 

we can expand (1 + /3mQ;~^)~^ according to 

00 

(1 + iima-^Y' = Y.{-\f^Va-^K (3.14) 

fe=0 

Upon exploiting the identity 

(i^)' = -\^^^ + ^^^^ + ^^-^^^ 

we find that W^"^ and W^^ do not contribute to C(0) by virtue of (2.50) and (2.51). The 
same holds for and , but the proof requires more intermediate steps, as follows. 

The term is given by 

-1 00 00 

Wlf = --(1 - 7)-' Y.i-lfP\k + 3)(A; + 5)(^ + 7) (m^ - m')m'^ a-^'-' . (3.16) 

Looking at even and odd values of /c, this suggests defining 

^ f m''^'a-^^'-\ G,^± m^'^^'a-''-', k = 0,l, 2..., (3.17) 

m=0 m=0 

G3 = E m^'^+^a-^'^-io, G4 = E m^'^+^a-^'^-^o. A; = 0, 1, 2... . (3.18) 

m=0 m=0 

Now Gi and 6*4 arc proportional to x"*^ and a;~^ respectively by virtue of (2.50), and hence 
do not contribute to C(0)- G2 contains weighted by a coefficient proportional to cos|, 
for all k, and hence does not contribute to C(0)- Last, G3 is such that its contribution 

Sgs{x; k) to C(0) reads 

1 _6 ^ r(5/2) r(4) r(5)\ 

5a3(a:;0) = --a; " 2^ + j = 0, (3-19) 
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jointly with infinitely many other relations along the same lines. 
The term W^^ is given by 

-1 oo oo 

w'^^ = T^(i - 7)-' E(-i)'/5'(^ + m + m + 5) E (^'^' - m'^')<^;;,'-'. (3.20) 

fe=0 m=0 

Here, too, we split the sum over k into sums over all even and odd values of k. We find 
therefore, exploiting (2.50) and (2.51), either contributions proportional to x~'^ and x~^, or 
terms weighted by cos |, or the contributions resulting from 

Hs^Y: m^'^'o^^'"-', ^ = 0, 1, 2, (3.21) 

m=0 

which occur with opposite signs for all k. 
The general formula for reads 

°° i'o '\2xdx) al 

oo / 1 I \ 3 



m 

r=0 



m=0 

4 2 

-2r- 

r=l r=0 



+ fe ^Y^j^^rrrf^J^ + ^ E i2,r-m^''a; 

r=l r=0 

= VF^^ + W'J + PF^^, (3.22) 



where jo,r,ji,r and j2,r are the coefficients in the polynomials flo, fli and respectively (see 
(2.37)-(2.39)), and negative powers of fe are expanded by exploiting 

(1 + fr = ti-^r^^^^f as / - 0. (3.23) 

Since 

1 2 00 

W^i^ = - o Eio,.(2r + 2)(2r + 4)(2r + 6) E (^"■^' - m'^^')a-'^-', (3.24) 

° r=0 m=0 

the basic formulae (2.50) and (2.51) imply a contribution to C(0) equal to 

^Ejo,. = 0. (3.25) 

Moreover, since 
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1 

OO 



r=l k=0 



k+r ^—k—r- 
m 



(3.26) 



m=0 



^^"^ = -^(l - 7)-' E^2,. E(-l)'(^ + + 2r + 2)(fc + 2r + 4){k + 2r + 6) 

° r=0 fe=0 



X ^(m 

m=0 



(3.27) 



repeated application of (2.50) and (2.51) yields contributions to C(0) equal to 



:i-7)-^(l + /5)-^E^V = 0, 



(3.28) 



r=l 



and 



1 



o(i - 7)-M E(-i)'(^ + E^v = 0, 



(3.29) 



\fe=0 / r=0 

respectively. The results (3.25), (3.28) and (3.29) are all vanishing because of the peculiar 
properties of the jo,r,ji,r ^-^d j2,r coefficients. 
Last, the general formula for reads 

= V (m^ - m) I — I \a~^ 
^ ^ ' \2xdx ] 



m=0 



E2r —2r 



r=0 



— r 
m 



r-=0 
3 



r=l 



m 



+ /e ^ E (^3,2rm 

r=0 



(3.30) 



where cro,2r, c"i,r, C2,r and (T3,2r are the coefficients in the polynomials u!o,u!i,u!2 and cus re- 
spectively (see (2.40)-(2.43)), and also fg^ is expanded by exploiting (3.23). Now the term 
W^^, which is the ^-independent part of (3.30), yields a non- vanishing contribution to C(0) 
equal to 



16 



while the terms W^^, W^'^ and W^^, resulting from ^j/^ ^ and fg ^ respectively, give 
vanishing contribution obtained as follows: 

^(l-7)-^(l + /?)-^(-Eai,,j ^OfromW^^^ (3.32) 



2 



^(1-7)-^ (E(-l)^ ^^^'f (- E '-3,2.) =0 from ly^^, (3.34) 

by repeated application of (2.50) and (2.51). 

By virtue of (3.4), (3.25), (3.28), (3.29), (3.31)-(3.34) we find 

for a massless Dirac field on the 4-ball, bearing in mind that also the eigenvalue condition 
(1.6) should be considered. Interestingly, such a ^(0) value in 4 dimensions is independent 
of 9, and agrees with the result in Ref. [12], where a massless spin-| field with half as many 
components as a Dirac field was instead considered. 

The proof of vanishing contributions to ^(0) from the infinite sums F^^G^ and can 
be made more systematic and elegant by remarking that a recursive scheme exists for which 

so that one only needs to look at F4 (see (3.10)), which can be evaluated exactly as a function 
of X for all k by exploiting the Euler-Maclaurin formula [3,12]. 
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IV. HEAT-KERNEL COEFFICIENTS IN GENERAL DIMENSION D 



Our aim in this section is to apply the formahsm in such a way that in principle all 
heat-kernel coefficients in any dimension d can be obtained. Therefore we will need the 
large-x expansion of (see (2.21) and (3.5)) 

oo A 



2x dx J 



2a, - {2u - 1) log(i/ + a,) + log(l + ^{t -1)) + Y1 



(4.1) 



^1 («-)^_ 

The first complication compared to d = 4 is that now, dealing with arbitrary dimension d, 
we need an arbitrary number of derivatives. This is easily generalized in some cases, i.e. 



m(2j-3)!! 1- 



2i 



a. 



2j 



in others at least in the form of a large- expansion. 

As we have seen in Sees. II and III, to deal with the Ap contributions, we need terms of the 

type 



j_d_ ^ ^ 

2idx) {1 + -f{t - iy)at ~ 2^ (1-7)' 

X 



V r 11" ^^^ + {u + t+p + 2j- 2)!! . 

' u\r{l) ^ (u + i + p-2)\\ ^ 



u=0 

The relevant case is j — 1 + d/2 and the contribution of each term to the zeta function is 
found by summing over n, taking the degeneracy into account. 

Let us now show how the general procedure works in the case of the ctj^-term. The 
contribution to ({1 + d/2,x'^) is 



^ ^ 2r 1 



d. 



- fd + n-2\ ( 1 dV^"\ 



+ d/2) ^0 V ^ 



^ {d-l)\\ - /rf + n - 2^ 
'2Wr(l + f)„^oV ^ y 



of which we need the large-x expansion. A simple expansion in inverse powers of x is not 
allowed; instead we employ a Mellin integral representation. As is clear from the above 
equation, the Barnes zeta function p5|-P^ 

Cb(s, a) = J2\^^ + (^y' = XI + + ••• + rnd-i)~% 

n=0 \ ^ / m=0 

will play a crucial role. We need to separate the u and x dependence in (4.2), more generally 
in expressions of the form 



CU, s) := E 



n 



n=0 



The V and x dependence is separated by employing for 3?c > 0, 

c+ioo 

e"^'* = / rfar(a)z/-2"t-°. 

27ri J 

c— ioo 

For 3fJs large enough we continue 



1 ^ (d + n-2\ , I 



c+ioo 



r(¥) 



n=0 



71=0 



n / 2712 

c— ioo 
c+ioo 

72 / 27[i 



C— 200 



The sum and integral may be interchanged upon choosing 3fJc > {j + d — l)/2 and we find 

c+ioo . / ^\ 



27rz 7 ' ' V 2 / ' \ 2 

c— ioo 



r(¥) 

On shifting the contour to the left we pick up the large-a; expansion of C{j,s). In order 
to find an/2, we are interested in the term that behaves as x"""^ and need to evaluate the 
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residue at a = (s + j — n — 2)/2. In all cases we encounter, the only relevant pole will come 
from and for these cases 



^^£=fcirri + ^^ RescJ.-n-2,^ 



2r(^) 



s) — 1 + - Res Cb s - n - 2, - x"""^ + irrelevant . (4.3) 



2/ " \ '2 



^From here, e.g., one obtains 



4v^ V 2 J '2j t{i + I) 

The procedure just outlined can be applied to all terms in (4.1). The following list sum- 
marizes for each term on the left the contribution to the heat-kernel coefficient an/2 on the 
right {6 and —6 contributions are summed): 

dg ( d — n — l\ ^ . ( . d\ 
(2a. - 2z.log(z. + a.)) ^ [^—^ j Res ^ (^ci - 1 - - j , 

1/ \ ds ^fd — n + l\ d\ 

log(. + a.)-^-^=^^r(^^^j ResCs[d-n,-j, 



log(l + 7(t - 1)) jT (^^^ (cosh-^- 9 - l) Res (s (d 



d\ 
' 2 / ' 



d^ ( d\ 

Res C,B\d + p — n, — \ x 



al{l + -i{t - l)y 4e'^cosh'^ ^ V '2 



r 



d+i+p—n 



r(^) 



^ + 1 / d + i+p-n 1 p + « ^ ,2 
3 2 1 2 ' 2' 2 2 2 



+ 6 ^-6. 

For the calculation of heat-kernel coefficients, note that only a finite number of terms con- 
tributes. The poles of the Barnes zeta function are located at s = 1, ...,d — l, and depending 
on the values of n and p only a finite number of terms needs to be evaluated. In general, to 
evaluate a„/2, we need to include all terms up to p = n — 1 [12]. 

The above results resemble very much the structure of the results found for different 
boundary conditions given in [^,^^-^. In particular, a reduction of the analysis from 



the ball to the sphere (in form of the Barnes zeta functions) has been achieved. Indeed, 
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instead of using the presented algorithm we could equally well have used the contour integral 
method developed in p3| , |29| , |30| . The starting point for the zeta function associated with the 
eigenvalues from (1.7) in this approach reads 

n=0 ^ 

the contour 7 enclosing all eigenvalues of (1.7). One then uses the uniform asymptotic 
expansion in order to extract the pieces that can contribute to the heat kernel coefficients. 
Performing the fc-integrals, results analogous to the above are found and final answers, of 
course, agree. 

Given the explicit results in the above list where all ingredients are known, the algo- 
rithm can be cast in a form suitable for application of Mathematica. As far as this process 
is concerned, some remarks are in order. We have presented the results in terms of hyper- 
geometric functions, and as far as we can see keeping rf, n arbitrary this is the best one can 
do. However, as soon as one considers particular values of d and n, the hypergeometric 
function 3F2 "collapses" to 2-^1, which, at the particular values needed, is simply given as 
an algebraic combination of hyperbolic functions. For example one has 

, 9 1 tanh 9 arcsin tanh 9 

2Fi(l,l,l/2,tanh^0) = ^ + — ^— — 

^ ' ' ^ ' ' l-tanh^^ (1 - tanh^ ^)3/2 

= cosh^ 9{1 + arcsin tanh 9 sinh 9). 

Mathematica will not always replace automatically the hypergeometric functions by this 
kind of hyperbolic combinations. Since this is essential for further simplifications of final 
answers, the implementation of some of the Gauss relations is necessary. We have used 

2Fi(« + l,/5+ 1,7+ l,z) = , {7 2Fi{a,(3,-i,z) - (7-«) 2Fi(«, /? + 1, 7 + 1, -2)} , 

7 2-Fi(a, 7, 2;) = (7 - a) 2Fi{a, /5, 7 + 1, z) + a 2-Fi(a + 1, /3, 7 + 1, ^). 

These relations guarantee that ultimately all hypergeometric functions are given in very 
explicit terms and that huge simplifications can be performed explicitly. In Sec. V we have 
summarized our findings in c? = 2, 4, 6 dimensions giving final results up to the coefficient 
ad/2- 
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V. LIST OF HEAT-KERNEL COEFFICIENTS 



We list hereafter the general results we have obtained. The lower coefficients have quite 
a simple form in all dimensions and the leading three coefficients are as follows: 

oo = — — — Tv> (5-1) 

(i + f ) 



a. - ^''-f: + U (l, h (tanh.)^l 

-(rf - 1)2^1(1, ^4 (tanh^f^j. (5.3) 

Moreover, to show the applicability of our algorithms to arbitrary dimensions and in principle 

to any coefficient we give the following collection of results. 

d=2: 

ao = ^, (5.4) 

ai/2 = ^(cosh^-l), (5.5) 

«i = (5.6) 
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d=4: 



ao = |, (5.7) 
ai/2 = ^((coshe)3-l), (5.8) 



ai = -^cosh2^, (5.9) 
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a3/2 = (sech (^\\ (15 + 20 cosh ^- 11 cosh 2^), 



(5.10) 



02 



114 

720 



(5.11) 



d=6: 



ao 



384' 



(5.12) 



0-1/2 



^((00*9)^-1 



(5,13) 



do 



ai 



384 



;-2 + 6cosh2^ + cosh4^), 



(5.14) 



ay 2 = fsech (^X) (153 + 212 cosh ^ + 35 cosh 2^ - 32 cosh 4^ - 8 cosh 50) , (5.15) 

98304 \ \ 2 / / 



0-2 = ^ cosh 26', 
9d 



(5.16) 
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05/2 = \ sech ( - 1 1 (311902 + 495474 cosh + 172792 cosh 20 

' 805306368 V \2// 

+ 14845 cosh 30 - 21590 cosh 40 - 2159 cosh 50) , 



(5.17) 



191c/, 



as 



(5.18) 



60480 

Of course, the result (5.11) for 02 in dimension four agrees with Eq. (3.35), upon bearing 
in mind that dg is then equal to 4. For = the results agree with the results found 



previously in [18,34,30 



VI. CONCLUDING REMARKS 



Motivated by quantum cosmology and the problems of quark confinement, we have stud- 
ied heat-kernel asymptotics for the squared Dirac operator on the Euclidean ball, with local 
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boundary conditions (1.1) leading to the eigenvalue conditions (1-6) and (1.7). We have first 
proved that on the 4-ball the ({0) value is 6'- independent. Furthermore, arbitrary values of d 
have been considered, and several explicit formulae for heat-kernel coefficients in dimension 
(i = 2, 4, 6 have been obtained in Sees. IV and V. Interestingly, ad/2 is always 6'-independent, 
while several other heat-kernel coefficients depend on 6 through hyperbolic functions and 
their integer powers. 

As far as we can see, the key task is now the analysis of heat-kernel asymptotics with 
local boundary conditions (1.1) on general Riemannian manifolds {M,g) with boundary 
dM. One has then to consider the smooth function / G C°°{M) mentioned after Eq. (1.13), 
which is replaced by 

a„/2(/, P, B) = c„/2(/, P) + 6„/2(/, P, B). (6.1) 

The interior part c„/2 vanishes for all odd values of n, whereas the boundary part only 
vanishes if = 0. The interior part is obtained by integrating over M a linear combination of 
local invariants of the appropriate dimension, where the coefficients of the linear combination 
are universal constants, independent of d. Moreover, the boundary part hn/2 is obtained upon 
integration over dM of another linear combination of local invariants. In that case, however, 
the structure group is 0{d ~ 1), and the coefficients of linear combination will depend on d 
and 6 [EH] and so they will be universal functions, as it happens if the boundary operator 



involves tangential derivatives |]3T| , |35| -p7[] . This is indeed the case for the boundary condition 
(1.1). To see this define 



and introduce the "projections" 



_ ■ 5 

X = ^ 7 7m 



n± = ^(i±x) 



In the bulk of our article we considered the operator P with domain 

domain(P) = E C°^{V) : U^i/j Ui © ^-D^P\g^^J = 0} 
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We calculate that 



n_DV \dM = (n;v^ + 7m7an*_v,) n+T/- 



dM 



a being a tangential index, which for Hermitean I\.± [9 = 0) would reduce to standard 
mixed boundary conditions. However, as is easily seen, this is not the case for 7^ and 
tangential derivatives occur in the boundary conditions such that the boundary conditions 
considered could be termed of mixed oblique type. It is thus expected, that the general 
form of an/2 contains all possible local invariants built from /, Riemann curvature R%cd of 
M, bundle curvature Qab (in case a gauge theory, with vector bundle over M, is studied), 
extrinsic curvature Kij of dM, endomorphism E (i.e. potential term) coming from the 
differential operator P, combinations of 7-matrices coming from the boundary operator, 
and the covariant derivatives of all these geometric objects, eventually integrating their 
linear combinations over M and dM [^p|,p|,|38| . All these local invariants are multiplied by 



universal functions which might depend on d and 6. As a next step, the presented special 
case calculation together with various other ingredients such as conformal transformations 
index theory [^], redefinition of the covariant derivative [^| will serve to find results 



valid for arbitrary Riemannian manifolds and bundle curvatures. 
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APPENDIX: 



The function T on the left-hand side of Eq. (1.7) is the product of the entire functions 
(i.e. functions analytic in the whole complex plane) 

^\ = Jn+d/2-l — Jn+d/2 and T2 = Jn+d/2-l + Jn+d/2-i 

which can be written in the form 

r^{k) = ^.k^+'^/^-'e^^^'^ ff fi - (Al) 



j=l V '^ij 

In Eqs. (Al) and (A2), 71 and 72 are constants, gi and g2 are entire functions, the fii are 
the zeros of J-'i and the z/j are the zeros of JF2. The general theory described in Ref. E2 



tells us that J-'i and JF2 are entire functions whose canonical product has genus 1. In other 
words, by virtue of the asymptotic behaviour of the eigenvalues, one finds that 

00 2 00 2 



00, 



=1 1/^ 



whereas and Y^^i jXTp are convergent. This is why the exponentials and 

must appear in Eqs. (Al) and (A2), which are called the canonical-product representations 
of JFi and JF2. The genus of the canonical product for JF^ is the minimum integer h such 
that converges, and similarly for F2, replacing /ij with z/j. If the genus is equal 

to 1, this ensures that no higher powers of ^ and ^ are needed in the argument of the 
exponential. Moreover, even for non-vanishing values of 6*, it remains true that the zeros of 
J^i are minus the zeros of J^2'- f^i = for all i [12]. Hence one finds eventually 

00 / p\ 

jr{k) = 7fc2(n+<i/2-i) (l-j^j, (A3) 

where 7 = 7172, /i^ are the positive zeros of J-'{k), and the sum {gi + g2){k) can be shown 
to vanish exactly as in Sec. IV of Ref. [12]. 
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In our paper we use uniform asymptotic expansions of regular Bessel functions and 
their first derivatives J^. On making the analytic continuation x ^ ix and then defining 
= + x^, one can write 

\/2tt 



Jlitx) ~ ML^ay2e°^e-'^^°^(^+"'')S2, (A5) 

V 277 



where the functions Si and S2 admit the asymptotic expansions 

00 00 
Si ~ ^Mfc(z//a^)/z/'', S2 ~ ^t;fc(z//a^)/ 

k=0 k=0 



valid uniformly in the order z/ as |x| cxd. The functions Uk and Vk are polynomials, given 



by Eqs. (9.3.9) and (9.3.13) on page 366 of Ref. |2| . 
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